Reversal-Field Memory in the Hysteresis of Spin Glasses 
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We report a novel singularity in the hysteresis of spin glasses, the reversal-field memory effect, 
which creates a non-analyticity in the magnetization curves at a particular point related to the 
history of the sample. The origin of the effect is due to the existence of a macroscopic number of 
"symmetric clusters" of spins associated with a local spin-reversal symmetry of the Hamiltonian. 
We use First Order Reversal Curve (FORC) diagrams to characterize the effect and compare to 
experimental results on thin magnetic films. We contrast our results on spin glasses to random 
magnets and show that the FORC technique is an effective "magnetic fingerprinting" tool. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 
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The non-equilibrium behavior of random magnets 
and spin glasses is an intensely studied field, pos- 
ing formidable theoretical and experimental challenges 
directly for magnetic systems, and also serving as 
paradigms for other fields. Concepts developed for ran- 
dom magnets such as glassy phases, droplet and replica 
theories, as well as aging have subsequently been applied 
to fields as diverse as structural biology, geology, and 
even financial analysis. 

The slow and complex time dependence of various cor- 
relators is a hallmark of such systems. Several aspects 
of this non-equilibrium dynamics have already been de- 
scribed in great detail for spin glasses Q. Hysteresis is 
one of the most central of these phenomena Q , yet while 
many basic features are qualitatively understood || [| , 
theoretical descriptions of hysteresis even in the simplest 
spin-glass models are in their early stages J(| . Hysteresis 
in magnetic systems has a host of practical applications 
including magnetic recording and sensors, but a less than 
complete understanding at a fundamental level 

In this paper, we present a detailed study of several 
new aspects of hysteresis in two of the most commonly 
studied models of disordered magnets. One is the Ran- 
dom Field Ising Model (RFIM) , which has been shown to 
describe successfully many of the relevant aspects of hys- 
teresis The second is the Edwards- Anderson Ising 
Spin Glass (EASG), which, unlike the RFIM, contains 
frustration, a phenomenon known to introduce a whole 
new level of complexity in disordered systems. Accord- 
ingly, we show that the hysteretic properties of the EASG 
can be significantly different from those of the RFIM. 

Our first important observation is of a novel memory 
effect in the hysteresis of the EASG that emerges when 
the magnetic field is first decreased from its saturation 
value and then increased again from some reversal field 
Hr. We find that the EASG exhibits a singularity at 
the negative of the reversal field, —Hr, in the form of 
a kink in the magnetization of the reversal curve. By 
calculating a suitable overlap function, we demonstrate 
that the microscopic origin of the effect is due to a macro- 



scopic number of "symmetric clusters" . In these clusters 
the central spins flip after all spins on the cluster surface 
have flipped. Therefore, the central spins experience an 
effective local field which is symmetric with respect to 
the change of direction of the external field. 

This reversal-field memory effect can be even more 
precisely characterized with the recently introduced 
First Order Reversal Curve (FORC) method @. As 
we shall demonstrate, the FORC technique provides a 
uniquely sensitive characterization of hysteretic systems 
and specifically of the difference between the hysteretic 
behavior of the RFIM and the EASG. The sharp kink of 
the minor loops of the EASG is captured as a profound 
horizontal ridge in FORC diagrams, indicative of a broad 
range of effective coercivities in the system, but a rather 
narrow range of biases. In contrast, despite exhibiting a 
major hysteresis loop rather similar to that of the EASG, 
the RFIM shows a strikingly different FORC diagram, 
characterized by a well developed vertical feature reflect- 
ing a rather narrow range of effective coercivities and a 
broad range of biases. 

Finally, we determine experimentally the reversal 
curves and FORC diagram of a magnetic thin film. Ex- 
perimentally, the reversal curves only show smoothed 
kinks around —Hr. However, the FORC diagram of the 
data reveals a profound horizontal ridge, signaling the 
presence of a reversal field memory in these films. This 
experimental result further highlights the usefulness of 
the FORC technique as a powerful method which cap- 
tures the detailed behavior of hysteretic systems. 

The Hamiltonian of the EASG is given by 



Ti — JijSiSj — H Si 



(1) 



(id) 



Here Si = ±1 are Ising spins on a square lattice of size 
N = LxLin two dimensions with periodic boundary con- 
ditions. The exchange couplings Jij arc random nearest- 
neighbor interactions chosen according to a Gaussian dis- 
tribution with zero mean and standard deviation unity, 
and H is the external magnetic field. We simulate the 
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FIG. 1: Reversal curve (solid line) and major hysteresis loop 
(dotted line) for a two-dimensional (2D) EASG with 10 4 spins 
and Hr — —2.28. In the inset a kink is seen around —Hr. In 
all figures the error bars are smaller than the symbols. 

zero temperature dynamics of the EASG by changing the 
external field H in small steps, first downward from posi- 
tive saturation and then upward from a reversal field Hr. 
After each field step, the effective local field hi of each 
spin Si is calculated: 

hi = J2 J a s 3- H - ( 2 ) 

3 

A spin is unstable if hiSi < 0. We then flip a randomly 
chosen unstable spin and update the local fields at neigh- 
boring sites and repeat this procedure until all spins are 
stable. 

Figure [j] (solid line) shows the average of 10 3 reversal 
curves, all with the same Hr, but different disorder real- 
izations. The area around —Hr is enlarged in the inset 
and shows a "kink." The presence of any such sharp fea- 
ture in a disordered system, especially of finite size and 
after disorder averaging, is quite remarkable. 

The change of slope at the kink can be characterized 
by measuring the slope of the magnetization curves to 
the left and right of —Hr, and comparing the difference 
A(dM/dH) with the average (dM/dH) ave (see Fig. |). 
The slope changes abruptly by as much as 30% as the 
field H = —Hr is passed, creating the kink. With our 
parameters the kink is present in the range of reversal- 
field values -4.0 < Hr < -1.5. 

In an effort to understand the microscopic origin of 
reversal-field memory, we first describe this effect within 
a phcnomenological approach to hysteretic systems, the 
Preisach model H . In the Preisach model a magnetic sys- 
tem is described as a collection of independent two-state 
(±1) switching units, or "hysterons". Unlike Ising spins, 
which always align with their local field, the hysteron's 
state changes from —1 to +1 at a field iff, + H c , different 
from the field Hi, — H c , required to switch the hysteron 
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FIG. 2: The difference (squares) and average (circles) of the 
left and right derivatives at — Hr for the EASG. 

from +1 to —1. Different systems are distinguished by 
their different distributions p(Hb,H c ) of hysterons of a 
given bias H^ and coercivity H c . Here p(Hi,,H c ) is the 
so-called "Preisach function" . 

An intuitive picture can be obtained by first consid- 
ering symmetric hysterons, having no bias, i.e. Hi, = 0. 
Starting from a fully "up" polarized state and decreasing 
the field to a negative Hr switches down all symmet- 
ric hysterons with H c < \Hr\. Reversing the direction of 
the sweep and increasing the field from Hr to — Hr along 
a reversal curve switches back every switched hysteron. 
Thus at H = —Hr saturation is reached, creating a kink 
in the magnetization. Symmetric hysterons therefore give 
rise to reversal-field memory. However, this memory ef- 
fect will be detectable only if the number of symmetric 
hysterons is macroscopic. This happens if p(Hb,H c ) has 
a Dirac delta singularity at Hi, = and H c = \Hr\. As 
the kink is observed in a range of Hr values, the singu- 
larities of the Preisach function form a horizontal ridge 
along the Hj, = axis for the corresponding range of 
H c = \Hr\ values. 

Next we move beyond phenomenological approaches, 
but keep the insight gained from the Preisach model. We 
carry over the concept of symmetric hysterons as "sym- 
metric clusters" of the strongly interacting spins of the 
EASG. A spin Si belongs to a symmetric cluster if Si 
flips down only after all its neighbors have flipped down, 
and during the reverse sweep Si flips up again only after 
all its neighbors have. Therefore, this central spin Si ex- 
periences an effective local field which is symmetric with 
respect to the change of direction of the external field, in 
analogy to a symmetric hysteron. 

Spins possessing local spin-reversal symmetry are can- 
didates for symmetric hysterons. By local spin-reversal 
symmetry we mean that the local field hi, felt by Si 
[Eq. ph], is perfectly reversed if the external field H is 
reversed and all spins coupled to Si are reversed as well. 
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FIG. 3: Overlap function q of the spins flipping at Hr and 
at H > Hr, for = —2.28. The insets show data for 
Hr — —0.40 and Hr = —5.60. The arrow in the inset marks 
H = 0.40. 

Every spin of the EASG has local spin-reversal symme- 
try. However, in a glassy system the spin configurations 
depend on the history of the sample. Therefore, at — Hr 
the neighbors of most spins do not necessarily point in a 
direction opposite of their direction at Hr, and thus most 
EASG spins do not belong to symmetric clusters. Hence 
the model Hamiltonian possessing a local spin-reversal 
symmetry is a necessary but not sufficient condition for 
having symmetric clusters. 

To see a macroscopic kink it has to be shown that the 
density of symmetric clusters is finite. A lower bound 
on their density is obtained by considering the simplest 
symmetric cluster: two strongly coupled spins, weakly 
coupled to their six neighbors in a 2D lattice. The switch- 
ing field of each spin is determined by these couplings. 
The outer spins will switch before the inner spins if their 
couplings are restricted by appropriate inequalities, con- 
fining the couplings to finite intervals. The density of 
symmetric clusters is obtained by integrating the prod- 
uct of the distributions of the couplings over these fi- 
nite intervals. With unbounded coupling distributions, 
e.g. Gaussian, the product of the distributions is finite 
over the finite integration intervals, thus the resulting 
density is finite as well. 

As a further evidence for the macroscopic number of 
symmetric clusters, we define an overlap function q be- 
tween the spins which flip at Hr and the spins which flip 
at H > H R : 

i 

[Si(H + 5) - Si{H)] . (3) 

Here 6 is the field step. In Fig. |3| we show the overlap 
q(H) for H R = -2.28. The large peak at H = +2.28 
indicates that a macroscopic number of spins which have 




FIG. 4: FORC Diagram of the EASG. Note the ridge along 
the H c axis. 

flipped at Hr, also flip at —Hr. This in turn means that 
there are a macroscopic number of symmetric clusters. 
The insets show a much smaller number of symmetric 
clusters at Hr = —0.40 and Hr = —5.60, values outside 
the peak of Fig. |[ 

To characterize reversal-field memory further, we 
adapt a new tool developed for analyzing experimen- 
tal data of hysteretic systems Q]. A family of First 
Order Reversal Curves (FORCs) with different Hr is 
generated, with M(H,Hr) denoting the resulting mag- 
netization as a function of the applied and reversal 
fields. Computing the mixed second order derivative 
p{H,Hr) = -(l/2)[d 2 M/dHdH R ] and changing vari- 
ables to H c = (H — H R )/2 and H b = (H + H R )/2, the 
local coercivity and bias, respectively, yields the "FORC 
distribution" p(Hf,,H c ). For phenomenological Preisach 
models, the FORC distribution is equal to the Preisach 
function. However, FORC distributions are more gen- 
eral, because they are extracted from numerical or ex- 
perimental data, and thus are model independent. 

Figure g shows the FORC diagram of the EASG. The 
ridge along the H c axis in the range 1.5 < H c < 4.0 
corresponds to the peak of Fig. |^, representing the kinks 
of Fig. |l|. Thus FORC diagrams capture the reversal-field 
memory effect in the form of a ridge along the H c axis. 

To demonstrate that local spin reversal symmetry of 
the Hamiltonian is necessary for reversal field memory 
to be present, we study the Random Field Ising Model 
(RFIM) H §. In this model J tj = 1 and the disor- 
der is introduced through random local fields chosen ac- 
cording to a Gaussian distribution with zero mean and 
standard deviation A. Direct inspection reveals that the 
RFIM does not possess a local spin-reversal symmetry. 
Therefore, the RFIM cannot have symmetric clusters and 
should not exhibit a reversal-field memory. This is con- 
firmed by our simulations: a typical RFIM reversal curve 
shown in the right panel of Fig. || has no kink at —Hr. 

Not only do the two models differ in the local spin 
symmetry, the EASG possesses frustration which might 
give rise to hysteretic phenomena that are qualitatively 
rather different than in the the RFIM. To explore this 
possibility, we show in the left panel of Fig. | the FORC 
diagram for the RFIM for a disorder A = 4.0. While the 
major hysteresis loop of the RFIM is very similar to that 
of the EASG, the FORC distribution is qualitatively dif- 
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FIG. 5: FORC Diagram of the 2D RFIM. Data for A = 4.0, 
TV = 10 4 spins and 10 3 disorder realizations (same scale as 
in Fig. /Q). Note the clear differences to the FORC diagram 
of the EASG (Fig. ^). The right panel shows the major hys- 
teresis loop (dotted line) and a reversal curve (solid line) for 
Hr — —1. The inset shows an enlarged view of the region 
around — Hr. 



ferent: it exhibits a predominantly vertical feature. The 
distribution of random fields of the RFIM introduces a 
large range of biases for the spins, with little variation 
in the local coercivity. This expectation is confirmed by 
simulations of the RFIM with various disorder distribu- 
tions, which show that the vertical cross section of the 
vertical feature mirrors the shape of the random field dis- 
tribution. 

Finally we demonstrate the existence of reversal-field 
memory in experimental systems. We study thin films of 
well-dispersed single-domain magnetic Co-7-Fe203 par- 
ticles provided by Kodak Inc. We determine both the 
individual reversal curves and the FORC diagram of the 
system. While the reversal field memory kinks in — Hr 
are somewhat smoothed, the FORC diagram clearly ex- 
hibits the horizontal ridge associated with the reversal- 
field memory effect (Fig. ^). This striking similarity be- 
tween the experimentally determined FORC diagram of 
the Co-7-Fe203 films and the numerically determined 
FORC diagram of the EASG indicates not only that Co- 
7-Fe203 films exhibit reversal-field memory but also that 
frustration may be a component of the physics of the Co- 
7~Fe203 films. 

In conclusion, we have reported a novel reversal-field 
memory effect in the EASG that manifests itself as a 
sharp kink in first order reversal curves and also as a 
sharp ridge on the zero bias axis of FORC diagrams. 
We suggest the microscopic origin of the effect is the 
presence of a macroscopic number of "symmetric clus- 
ters," and prove this by computing a suitable overlap 




FIG. 6: Experimental FORC diagram of a Kodak sample. 
Note the similarity to the FORC diagram of the EASG shown 
in Fig. | 



function. We further show that reversal field memory is 
absent from the RFIM, which does not exhibit symmet- 
ric clusters. While the hysteresis loops of the EASG and 
the RFIM are remarkably similar for corresponding pa- 
rameters, their FORC diagrams are profoundly different, 
establishing that the FORC method is a powerful diag- 
nostic tool for capturing the sensitive details of hysteretic 
systems such as spin glasses and random magnets. The 
FORC diagrams of several magnetic thin films exhibit 
a profound ridge indicative of the reversal-field memory 
effect in experimental systems. Simulations on more real- 
istic magnetic models which include dipolar interactions 
have also shown the reversal-field memory effect. This 
suggests that the reversal-field memory is not specific to 
the EASG, but it a robust result for a large class of the- 
oretical models and experimental systems. 
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